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Abstract 

In this paper, applying a canonical system with field rotation parameters and using 
geometric properties of the spirals filling the interior and exterior domains of limit 
cycles, we solve first the problem on the maximum number of limit cycles surround- 
ing a unique singular point for an arbitrary polynomial system. Then, by means 
of the same bifurcationally geometric approach, we solve the limit cycle problem 
for a general Lienard polynomial system with an arbitrary (but finite) number of 
singular points. This is related to the solution of Hilbert's sixteenth problem on 
the maximum number and relative position of limit cycles for planar polynomial 
dynamical systems. 
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1 Introduction 



We will consider Lienard equations of the form 

x + f(x) x + g(x) = 0. (1.1) 

In the phase plane, the representation of equation (1.1) is given by the dy- 
namical system 

x = y, y = -g(x) - f(x)y. (1.2) 
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There are many examples in the natural sciences and technology in which this 
and related systems are applied [1] [13] . Such systems are often used to model 
either mechanical or electrical, or biomedical systems, and in the literature, 
many systems are transformed into Lienard type to aid in the investigations. 
They can be used, e. g., in certain mechanical systems, where f(x) represents a 
coefficient of the damping force and g(x) represents the restoring force or stiff- 
ness, when modeling wind rock phenomena and surge in jet engines [2], [TO] . 
Such systems can be also used to model resistor-inductor-capacitor circuits 
with non-linear circuit elements. Recently, e.g., the Lienard system (1.2) has 
been shown to describe the operation of an optoelectronics circuit that uses a 
resonant tunnelling diode to drive a laser diode to make an optoelectronic vol- 
tage controlled oscillator [12]. There are also some examples of using Lienard 
type systems in ecology and epidemiology [9]. 

In this paper, we suppose that system (1.2), where f(x) and g(x) are arbitrary 
polynomials of x, has an anti-saddle (a node or a focus, or a center) at the 
origin and write it in the form 

x = y, y — — x (1 + /3i x + . . . + /3 2 i x 21 ) +y(a + a 1 x + ... + a 2 k x 2k ). (1.3) 

Note that for g(x) = x, by the change of variables X = x and Y = y + F(x), 
where F(x) = Jq f{s) ds, (1.3) is reduced to an equivalent system 

X = Y-F(X), Y = -X (1.4) 

which can be written in the form 

i = V, y = ~x + F(y) (1.5) 

or 

x = y, y = -x + j 1 y + -f 2 y 2 + i3y 3 + --- + i2ky 2k + i2k+iy 2k+1 - (1.6) 

Therefore, we can conclude that our previous results [H], [15] agree with the 
conjecture of [7] on the maximum number of limit cycles for the classical 
Lienard polynomial system (1.6). In [TJ], |15j . we have presented a solution of 
Smale's thirteenth problem JT3] proving that the Lienard system (1.6) with a 
polynomial of degree 2k + 1 can have at most k limit cycles. In [16j-[2T], we 
have also presented a solution of Hilbert's sixteenth problem in the quadratic 
case of polynomial systems proving that for quadratic systems four is really the 
maximum number of limit cycles and (3: 1) is their only possible distribution. 
We have established some preliminary results on generalizing our ideas and 
methods to special cubic, quartic and other polynomial dynamical systems 
as well. In [22], e.g., we have constructed a canonical cubic dynamical system 
of Kukles type and have carried out the global qualitative analysis of its spe- 
cial case corresponding to a generalized Lienard equation. In particular, it has 
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been shown that the foci of such a Lienard system can be at most of second 
order and that such system can have at most three limit cycles in the whole 
phase plane. Moreover, unlike all previous works on the Kukles-type systems, 
global bifurcations of limit and separatrix cycles using arbitrary (including as 
large as possible) field rotation parameters of the canonical system have been 
studied. As a result, a classification of all possible types of separatrix cycles 
for the generalized Lienard system has been obtained and all possible distri- 
butions of its limit cycles have been found. In [23], [24], we have completed the 
global qualitative analysis of a planar Lienard-type dynamical system with a 
piecewise linear function containing an arbitrary number of dropping sections 
and approximating an arbitrary polynomial function. In [25] [29], we have also 
completed the global qualitative analysis of a quartic dynamical system which 
models the dynamics of the populations of predators and their prey in a given 
ecological system. In [30] , we have studied of multiple limit cycle bifurcations 
in the well-known FitzHgh-Nagumo neuronal model. 

We use the obtained results and develop our methods for studying limit cycle 
bifurcations of polynomial dynamical systems in this paper as well. In Sec- 
tion 2, applying a canonical system with field rotation parameters and using 
geometric properties of the spirals filling the interior and exterior domains of 
limit cycles, we solve first the problem on the maximum number of limit cycles 
surrounding a unique singular point for an arbitrary polynomial system. Then, 
in Section 3, by means of the same bifurcationally geometric approach, we solve 
the limit cycle problem for the general Lienard polynomial system (1.3) with 
an arbitrary (but finite) number of singular points proving that system (1.3) 
can have at most k + I limit cycles, k surrounding the origin and / surrounding 
one by one the other singularities of (1.3). This is related to the solution of 
Hilbert's sixteenth problem on the maximum number and relative position of 
limit cycles for planar polynomial dynamical systems. 



2 Limit cycles surrounding a unique singular point 

Consider first an arbitrary dynamical system 

x = P n (x,y,n 1} . . . ,Hk), y = Qn(x,y,fj, 1 ,...,fi k ), (2.1) 

where P n and Q n are polynomials in the real variables x, y and not greater 
than n-th degree containing k field rotation parameters, fj,i, . . . , and having 
an anti-saddle at the origin. We prove the following theorem. 

Theorem 2.1. The polynomial system (2.1) containing k field rotation pa- 
rameters and having a singular point of center type at the origin for the zero 
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values of these parameters can have at most k — 1 limit cycles surrounding the 
origin. 

Proof. Let all the parameters of (2.1) vanish and suppose that the obtained 
system 

x = P n (x, y, 0, . . . , 0), y= Q n {x, y, 0, . . . , 0) (2.2) 
has a singular point of center type at the origin. 

Input successively the field rotation parameters, /xi, . . . , /^., into this system 
(see [H], [IS]). Suppose, e.g., that > and that the vector field of the 
system 

x = P n (x,y,/xi,0, . . . ,0), y = Q n (x,y, fjL 1} 0, . . . , 0) (2.3) 

is rotated counterclockwise turning the origin into a stable focus under in- 
creasing the parameter /j,x [I], [TJ], [T5] . 

Fix yUi and input the parameter yu 2 into (2.3) changing it so that the field of 
the system 

x = P n (x,y, ^1,^2,0, ... ,0), y = Q n (x,y, fa, /j.2,0, . . . , 0) (2.4) 

would be rotated in the opposite direction (clockwise). Suppose this occurs 
for /i 2 < 0. Then, for some value of this parameter, a limit cycle will appear in 
system (2.4). There are three logical possibilities for such a bifurcation: 1) the 
limit cycle appears from the focus at the origin; 2) it can also appear from 
some separatrix cycle surrounding the origin; 3) the limit cycle appears from 
a so-called "trajectory concentration" . In the last case, the limit cycle is semi- 
stable and, under further decreasing it splits into two limit cycles (stable 
and unstable) one of which then disappears at (or tends to) the origin and the 
other disappears on (or tends to) some separatrix cycle surrounding this point. 
But since the stability character of both a singular point and a separatrix 
cycle is quite easily controlled [18] . this logical possibility can be excluded. 
Let us choose one of the two other possibilities: e.g., the first one, the so- 
called Andronov-Hopf bifurcation. Suppose that, for some value of /z 2 , the 
focus at the origin becomes non-rough, changes the character of its stability 
and generates a stable limit cycle, T\. 

Under further decreasing /z 2 , three new logical possibilities can arise: 1) the 
limit cycle Li disappears on some separatrix cycle surrounding the origin; 
2) a separatrix cycle can be formed earlier than Y\ disappears on it and then 
it generates one more (unstable) limit cycle, T 2 , which joins with Y\ forming 
a semi-stable limit cycle, r 12 , disappearing in a "trajectory concentration" 
under further decreasing /i 2 ; 3) in the domain D\ outside the cycle Li or in 
the domain Z} 2 inside Yx, a semi-stable limit cycle appears from a "trajectory 
concentration" and then splits into two limit cycles (logically, the appearance 
of such semi-stable limit cycles can be repeated). 
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Let us consider the third case. It is clear that, under decreasing fi 2 , a semi- 
stable limit cycle cannot appear in the domain D 2 , since the focus spirals filling 
this domain will untwist and the distance between their coils will increase 
because of the vector field rotation. By contradiction, we can prove that a 
semi-stable limit cycle cannot appear in the domain D\. 

Suppose it appears in this domain for some values of the parameters ^\ > 
and \i\ < 0. Return to system (2.2) and change the inputting order for the 
field rotation parameters. Input first the parameter /i 2 < : 

x = P n (x,y,/j, 2 ,0, ...,0), y = Q n (x,y,/j, 2 ,0,...,0). (2.5) 

Fix it under fj, 2 = n 2 - The vector field of (2.5) is rotated clockwise and the 
origin turns into an unstable focus. Inputting the parameter \i\ > into (2.5), 
we get again system (2.4) the vector field of which is rotated counterclockwise. 
Under this rotation, a stable limit cycle, I\, will appear from some separatrix 
cycle. The limit cycle I\ will contract, the outside spirals winding onto this 
cycle will untwist and the distance between their coils will increase under 
increasing Hi to the value fi*. It follows that there are no values of \i\ < and 
/4 > for which a semi-stable limit cycle could appear in the domain D\. 

The second logical possibility can be excluded by controlling the stability 
character of the separatrix cycle [IS]. Thus, only the first possibility is valid, 
i.e., system (2.4) has at most one limit cycle. 

Let system (2.4) have the unique limit cycle T\. Fix the parameters /ii > 0, 
H 2 < and input the third parameter, /i 3 > 0, into this system supposing that 
fi 3 rotates its vector field counterclockwise: 

x = P n (x,y,n 1} (i2,na } 0,...,0), y = Q n {x, y, pi, /x 2 , fi 3 , 0, . . . , 0). (2.6) 

Here we can have two basic possibilities: 1) the limit cycle Fl disappears at the 
origin; 2) the second (unstable) limit cycle, T 2 , appears from the origin and, 
under further increasing the parameter /i 3 , the cycle T 2 joins with I\ forming a 
semi-stable limit cycle, Ti2, which disappears in a "trajectory concentration" 
surrounding the origin. Besides, we can also suggest that: 3) in the domain D 2 
bounded by the origin and Ft, a semi-stable limit cycle, T23, appears from a 
"trajectory concentration" , splits into an unstable cycle, T 2 , and a stable cycle, 
r 3 , and then the cycles r 1; T 2 disappear through a semi-stable limit cycle, r 12 , 
and the cycle T 3 disappears through an Andronov-Hopf bifurcation; 4) a semi- 
stable limit cycle, T34, appears in the domain D 2 bounded by the cycles Ft, 
T 2 and, for some set of values of the parameters, /i*, [i 2 , fi 3 , system (2.6) has 
at least four limit cycles. 

Let us consider the last, fourth case. It is clear that a semi-stable limit cycle 
cannot appear either in the domain Di bounded on the inside by the cycle Ft 
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or in the domain D 3 bounded by the origin and T 2 because of the increasing 
distance between the spiral coils filling these domains under increasing the 
parameter yu 3 . To prove the impossibility of the appearance of a semi-stable 
limit cycle in the domain D 2 , suppose the contrary, i. e., for some set of values 
of the parameters, /i{ > 0, /x 2 < 0, and // 3 > 0, such a semi-stable cycle exists. 
Return to system (2.2) again and input first the parameters ^3 > 0, [i\ > 0: 

x = P n (x,y,fi 1 ,fi 3 ,0,...,0), y = Q n (x,y,ni,fi 3 ,0,...,0). (2.7) 

Fix these parameters under // 3 = /i 3 , /ii = fi\ and input the parameter /x 2 < 
into (2.7) getting again system (2.6). Since, by our assumption, this system 
has two limit cycles for /i 2 > 1^2, there exists some value of the parameter, 
/4 2 (a*2 < /4 2 < 0)i f° r which a semi-stable limit cycle, r 12 , appears in system 
(2.6) and then splits into a stable cycle, r 1; and an unstable cycle, T 2 , under 
further decreasing /i 2 . The formed domain D2 bounded by the limit cycles 
Ft, T 2 and filled by the spirals will enlarge, since, by the properties of a field 
rotation parameter, the interior unstable limit cycle T 2 will contract and the 
exterior stable limit cycle r\ will expand under decreasing /z 2 . The distance 
between the spirals of the domain D 2 will naturally increase, what will prevent 
the appearance of a semi-stable limit cycle in this domain for /i 2 < a4 2 - 

Thus, there are no such values of the parameters, /i* > 0, /i 2 < 0, /i 3 > 0, for 
which system (2.6) would have an additional semi-stable limit cycle. Therefore, 
the fourth case cannot be realized. The third case is considered absolutely 
similarly. It follows from the first two cases that system (2.6) can have at 
most two limit cycles. 

Suppose that system (2.6) has two limit cycles, I\ and T 2 , fix the parameters 
/ii > 0, fj,2 < 0, fj,3 > and input the fourth parameter, // 4 < 0, into this 
system supposing that /i 4 rotates its vector field clockwise: 

x = P n (x,y,Hi,...,nt,0,...,0), y = Q n (a;,y,/ii,...,/i4,0,...,0). (2.8) 

The most interesting logical possibility here is that when the third (stable) 
limit cycle, T 3 , appears from the origin and then, under preservation of the 
cycles r\ and T 2 , in the domain D 3 bounded on the inside by the cycle T 3 
and on the outside by the cycle T 2 , a semi-stable limit cycle, r 45 , appears 
and then splits into a stable cycle, r 4 , and an unstable cycle, T^, i.e., when 
system (2.8) for some set of values of the parameters, /i*, /i 2 , /Xg, £i 4 , has at least 
five limit cycles. Logically, such a semi-stable limit cycle could also appear in 
the domain D 1 bounded on the inside by the cycle r\ , since, under decreasing 
/x 4 , the spirals of the trajectories of (2.8) will twist and the distance between 
their coils will decrease. On the other hand, in the domain _D 2 bounded on 
the inside by the cycle T 2 and on the outside by the cycle F and also in the 
domain D4 bounded by the origin and T 3 , a semi-stable limit cycle cannot 
appear, since, under decreasing ^ 4 , the spirals will untwist and the distance 
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between their coils will increase. To prove the impossibility of the appearance 
of a semi-stable limit cycle in the domains D 3 and Di, suppose the contrary, 
i.e., for some set of values of the parameters, Hi > 0, \x\ < 0, nt > 0, and 
H% < 0, such a semi-stable cycle exists. 

Return to system (2.2) again, input first the parameters Hi < 0, Hi < and 
then the parameter Hi > : 

x = P n (x,y,Hi,ti 2 ,tJ.A,0, ■ ■ ■ ,0), y = Q n (x,y, pi, 112,114,0,..., 0). (2.9) 

Fix the parameters H4, //2 under the values n%, /4, respectively. Under increas- 
ing Hi, a separatrix cycle is formed around the origin generating a stable limit 
cycle, Ti. Fix Hi under the value Hi an d input the parameter /i 3 > into (2.9) 
getting system (2.8). 

Since, by our assumption, system (2.8) has three limit cycles for /i 3 < H31 there 
exists some value of the parameter hI 3 (0 < H3 3 < ^3) 101 which a semi-stable 
limit cycle, T 2 3, appears in this system and then splits into an unstable cycle, 
T2, and a stable cycle, T3, under further increasing ^3. The formed domain D 3 
bounded by the limit cycles r 2 , T 3 and also the domain Di bounded on the 
inside by the limit cycle I\ will enlarge and the spirals filling these domains 
will untwist excluding a possibility of the appearance of a semi-stable limit 
cycle there. 

All other combinations of the parameters Hi, Hit an d Hi are considered in 
a similar way. It follows that system (2.8) has at most three limit cycles. If we 
continue the procedure of successive inputting the field rotation parameters, 
H5, He, ■ ■ ■ , Hki into system (2.2), it is possible to conclude that system (2.1) 
can have at most k — 1 limit cycles surrounding the origin. The theorem is 
proved. □ 



3 Limit cycles of the general Lienard polynomial system 

By means of the same bifurcationally geometric approach, we will consider 
now the Lienard polynomial system (1.3). The study of singular points of 
system (1.3) will use two index theorems by H. Poincare, see pp. The definition 
of the Poincare index is the following pQ. 

Definition 3.1. Let S be a simple closed curve in the phase plane not passing 
through a singular point of the system 

x = P(x,y), y = Q(x,y), (3.1) 

where P(x, y) and Q(x, y) are continuous functions (for example, polynomials), 
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and M be some point on S. If the point M goes around the curve S in the 
positive direction (counterclockwise) one time, then the vector coinciding with 
the direction of a tangent to the trajectory passing through the point M is 
rotated through the angle 2tvJ (j = 0, ±1, ±2, . . .). The integer j is called the 
Poincare index of the closed curve S relative to the vector field of system (3.1) 
and has the expression 

. 1 rPdQ-QdP 



According to this definition, the index of a node or a focus, or a center is equal 
to +1 and the index of a saddle is —1. 

Theorem 3.1 (First Poincare Index Theorem). If N, Nf, N c , and C are 

respectively the number of nodes, foci, centers, and saddles in a finite part of 
the phase plane and N' and C are the number of nodes and saddles at infinity, 
then it is valid the formula 



N + N f + N c + N' = C + C' + l. (3.3) 



Theorem 3.2 (Second Poincare Index Theorem). If all singular points 
are simple, then along an isocline without multiple points lying in a Poincare 
hemisphere which is obtained by a stereographic projection of the phase plane, 
the singular points are distributed so that a saddle is followed by a node or a 
focus, or a center and vice versa. If two points are separated by the equator 
of the Poincare sphere, then a saddle will be followed by a saddle again and a 
node or a focus, or a center will be followed by a node or a focus, or a center. 

Consider system (1.3). Its finite singularities are determined by the algebraic 
system 

x{l + p lX + ... + (3 2l x 21 ) = 0, y = 0. (3.4) 

It always has an anti-saddle at the origin and, in general, can have at most 
21 + 1 finite singularities which lie on the x-axis and, according to Theorem 3.2, 
are distributed so that a saddle (or saddle-node) is followed by a node or a 
focus, or a center and vice versa. At infinity, system (1.3) has two singular 
points: a node at the "ends" of the x-axis and a saddle at the "ends" of the 
y-axis. For studying the infinite singularities, the methods applied in [1] for 
Rayleigh's and van der Pol's equations and also Erugin's two-isocline method 
developed in [TS] can be used (see [H], [To]). 

Following [18] . we will study limit cycle bifurcations of (1.3) by means of 
a canonical system containing field rotation parameters of (1.3) [T], [T5] . 
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Theorem 3.3. The Lienard polynomial system (1.3) with limit cycles can be 
reduced to the canonical form 

x = y = P(x,y), 
y = - x (1 + ft x ± x 2 + . . . + p 2 i-i x 21 ' 1 ± x 21 ) (3-5) 
+ y (a + x + a 2 x 2 + . . . + x 2fc_1 + a 2k x 2k ) = Q(x, y), 

where ft, ft, ... , P21-1 ar e fixed and ao, a 2 , ■ ■ ■ , a2fe ore /ie/<i rotation para- 
meters of (3.5). 

Proof. Let all the parameters cej, i — 0, 1, . . . , 2&, vanish in system (3.5), 

i; = y, y= -x (1 + ft x + ft x 2 + . . . + ft/ x 2i ), (3.6) 
and consider the corresponding equation 

dy -x (1 + ft x + ft x 2 + . . . + ft/ x 2 ') 



dx y 



F(x,y). (3.7) 



Since F(x, — y) = —F(x, y), the direction field of (3.7) (and the vector field of 
(3.6) as well) is symmetric with respect to the x-axis. It follows that for arbi- 
trary values of the parameters ft, j — 1, 2, . . . , 21, system (3.6) has centers as 
anti-saddles and cannot have limit cycles surrounding these points. Therefore, 
without loss of generality, all the even parameters ft of system (1.3) can be 
supposed to be equal, e. g., to ±1: ft = ft = ft = . . . = ±1. 

Let now all the parameters with even indexes and ft with odd indexes 
vanish in system (3.5), 

x — y, y = -x (1± x 2 ±. . .± x 2l )+y (a x x + a 3 x 3 + . . . + a 2 k-i x 2 ^ 1 ), (3.8) 

and consider the corresponding equation 

dy — x (1 ± x 2 ± . . . ± x 21 ) + y (an x + a 3 x 3 + . . . + a 2 k-i x 2fc_1 ) 
dx y (3.9) 

= G(x,y). 

Since G(—x,y) = —G{x,y), the direction field of (3.9) (and the vector field 
of (3.8) as well) is symmetric with respect to the y-axis. It follows that for 
arbitrary values of the parameters «i, 0:3, . . . , a 2 k-i system (3.6) has centers as 
anti-saddles and cannot have limit cycles surrounding these points. Therefore, 
without loss of generality, all the odd parameters 0;^ of system (1.3) can be 
supposed to be equal, e.g., to 1: a± — a 3 — . . . — a 2 k-i = 1- 
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Inputting the odd parameters fa, (3s, ... , fai-i into system (3.8), 



x = y = R(x,y), 

y = -x {1 + fax ± x 2 + fax 3 ± x 4 + . . . + (5 2 i-i x 21 - 1 ± x 21 ) (3.10) 
+ y(x + x 3 + ... + x 2fc " 1 ) = S(x, y), 

and calculating the determinants 

A ft = RS'h - SR'^ = -x 2 y, 
A fe = RS^ - SR' Pz = -x 4 y, 

Aa,., = RS'^ - SR!^ = -x 2l y, 

we can see that the vector field of (3.10) is rotated symmetrically (in opposite 
directions) with respect to the x-axis and that the finite singularities (centers 
and saddles) of (3.10) moving along the (except the center at the origin) 

do not change their type or join in saddle-nodes. Therefore, we can fix the odd 
parameters fa, fa, ... , fai-i in system (3.5), fixing the position of its finite 
singularities on the x-axis. 

To prove that the even parameters ao, ... , «2fc rotate the vector field 
of (3.5), let us calculate the following determinants: 

Aa = PQ' ao - QK Q = y 2 >o, 
a« 2 = PQ' a2 - QK 2 = * 2 y 2 > o, 



A« 2fc = PQ' a2k - QPL 2k = * 2k y 2 > o. 

By definition of a field rotation parameter PQ, [18], for increasing each of 
the parameters ao, a 2 , ... , a 2 k, under the fixed others, the vector field of 
system (3.5) is rotated in the positive direction (counterclockwise) in the whole 
phase plane; and, conversely, for decreasing each of these parameters, the 
vector field of (3.5) is rotated in the negative direction (clockwise). 

Thus, for studying limit cycle bifurcations of (1.3), it is sufficient to consider 
the canonical system (3.5) containing only its even parameters ao, a 2 , . . . , a 2 k 
which rotate the vector field of (3.5) under the fixed others. The theorem is 
proved. □ 

By means of the canonical system (3.5), let us study global limit cycle bifur- 
cations of (1.3) and prove the following theorem. 
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Theorem 3.4. The general Lienard polynomial system (1.3) can have at most 
k + I limit cycles, k surrounding the origin and I surrounding one by one the 
other singularities of (1.3). 

Proof. According to Theorem 3.3, for the study of limit cycle bifurcations 
of system (1.3), it is sufficient to consider the canonical system (3.5) contain- 
ing the field rotation parameters ao, a 2 , • • • , ct2k of (1.3) under the fixed its 
parameters ft, ft, ... , ft/_i. 

Let all these parameters vanish: 

± = y, y = -x(l± x 2 ±...± x 2l ) + y(x + x 3 + ... + x 2k ' 1 ). (3.11) 

System (3.11) is symmetric with respect to the y-axis and has centers as anti- 
saddles. Its center domains are bounded by either separatrix loops or digons of 
the saddles of (3.11) lying on the x-axis. If to input the parameters ft, ft, ... , 
ft/_i into (3.11) successively, we will get again system (3.10) the vector field 
of which is rotated symmetrically (in opposite directions) with respect to the 
x-axis. The finite singularities (centers and saddles) of (3.10) moving along 
the x-axis (except the center at the origin) do not change their type or join in 
saddle-nodes and the center domains will be bounded by separatrix loops of 
the saddles (or saddle-nodes) of (3.10) [T], [T8] . 

Let us input successively the field rotation parameters ao, a 2 , ■ ■ ■ , a 2 k into 
system (3.10) beginning with the parameters at the highest degrees of x and 
alternating with their signs (see [H], [IS]). So, begin with the parameter a 2 k 
and let, for definiteness, a 2 k > 0: 

x = y, 

y = -x (1 + ft X ± X 2 + ft X 3 ± X 4 + . . . + ft;_! X 2 '- 1 ± X 21 ) (3.12) 

+ y (x + x 3 + . . . + x 2fc ~ 1 + a 2 k x 2k ) . 

In this case, the vector field of (3.12) is rotated in the positive direction (coun- 
terclockwise) turning the center at the origin into a nonrough (weak) unstable 
focus. All the other centers become rough unstable foci [lj, [18J. 

Fix a 2 fc and input the parameter a 2 k~ 2 < into (3.12): 

x = y, 

y = -x (1 + ft X ± X 2 + ft X 3 ± X 4 + . . . + ft;_! X 2 ^ 1 ± X 2 ') (3.13) 

+ y (x + x 3 + . . . + a 2 k-2X 2h ~ 2 + x 2k ~ l + a 2 k x 2k ). 
Then the vector field of (3.13) is rotated in the opposite direction (clockwise) 
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and the focus at the origin immediately changes the character of its stability 
(since its degree of nonroughness decreases and the sign of the field rotation 
parameter at the lower degree of x changes) generating a stable limit cycle. All 
the other foci will also generate stable limit cycles for some values of a 2 k- 2 after 
changing the character of their stability Under further decreasing a 2 k- 2 , all 
the limit cycles will expand disappearing on separatrix cycles of (3.13) [I], [T8] . 

Denote the limit cycle surrounding the origin by Fx , the domain outside the cy- 
cle by Di, the domain inside the cycle by D 2 and consider logical possibilities 
of the appearance of other (semi-stable) limit cycles from a "trajectory con- 
centration" surrounding this singular point. It is clear that, under decreasing 
the parameter a 2 k- 2 , a semi-stable limit cycle cannot appear in the domain 
D 2 , since the focus spirals filling this domain will untwist and the distance 
between their coils will increase because of the vector field rotation [TJ], [T5] . 

By contradiction, we can also prove that a semi-stable limit cycle cannot 
appear in the domain D\. Suppose it appears in this domain for some values 
of the parameters a 2k > and a 2k _ 2 < 0- Return to system (3.10) and change 
the inputting order for the field rotation parameters. Input first the parameter 
«2fc-2 < : 

x = y, 

y = -x (1 + ft x ± x 2 + p 3 x 3 ± x 4 + . . . + P21-1 x 2 ^ 1 ± x 21 ) (3.14) 
+ y(x + x 3 + ... + a 2k _ 2 x 2k ~ 2 + x 2fc_1 ). 

Fix it under a 2 k-2 = a 2k-2- The vector field of (3.14) is rotated clockwise and 
the origin turns into a nonrough stable focus. Inputting the parameter a 2 k > 
into (3.14), we get again system (3.13) the vector field of which is rotated 
counterclockwise. Under this rotation, a stable limit cycle Fi will appear from 
a separatrix cycle for some value of a 2 k- This cycle will contract, the outside 
spirals winding onto the cycle will untwist and the distance between their coils 
will increase under increasing a 2 k to the value a 2k . It follows that there are 
no values of u 2k _ 2 < and a 2k > for which a semi-stable limit cycle could 
appear in the domain D\. 

This contradiction proves the uniqueness of a limit cycle surrounding the origin 
in system (3.13) for any values of the parameters a 2 k- 2 and a 2 k of different 
signs. Obviously, if these parameters have the same sign, system (3.13) has 
no limit cycles surrounding the origin at all. On the same reason, this system 
cannot have more than / limit cycles surrounding the other singularities (foci 
or nodes) of (3.13) one by one. 

Let system (3.13) have the unique limit cycle Ft surrounding the origin and 
I limit cycles surrounding the other antisaddles of (3.13). Fix the parameters 



12 



«2fe > 0, «2fc-2 < and input the third parameter, a2k-A > 0, into this system: 

x = y, 

y = - x (1 + ft x ± x 2 + ft x 3 ± x 4 + . . . + p 2 i-i x 21 - 1 ± x 21 ) (3.15) 

+ y (x + x 3 + . . . + a 2k -iX 2k ~ 4 + a 2k - 2 x 2k ~ 2 + x 2h ~ l + a 2k x 2k ). 

The vector field of (3.15) is rotated counterclockwise, the focus at the origin 
changes the character of its stability and the second (unstable) limit cycle, T 2 , 
immediately appears from this point. The limit cycles surrounding the other 
singularities of (3.15) can only disappear in the corresponding foci (because of 
their roughness) under increasing the parameter 0:2*1-4. Under further increas- 
ing a 2 k-4, the limit cycle T 2 will join with T\ forming a semi-stable limit cycle, 
1x2, which will disappear in a "trajectory concentration" surrounding the ori- 
gin. Can another semi-stable limit cycle appear around the origin in addition 
to Ti2? It is clear that such a limit cycle cannot appear either in the domain 
D\ bounded on the inside by the cycle Y\ or in the domain D3 bounded by 
the origin and Y 2 because of the increasing distance between the spiral coils 
filling these domains under increasing the parameter a 2 k-i [Ej, [15] . 

To prove the impossibility of the appearance of a semi-stable limit cycle in the 
domain D 2 bounded by the cycles Y\ and Y 2 (before their joining), suppose the 
contrary, i. e., that for some set of values of the parameters, a\ k > 0, a 2k _ 2 < 0, 
and Q!2A;-4 > O5 suc h a semi-stable cycle exists. Return to system (3.10) again 
and input first the parameters a 2 k-4 > and a 2 k > : 

x = y, 

y = -x (1 + ft X ± X 2 + P 3 X 3 ± X 4 + . . . + 021-1 X 21 - 1 ± X 21 ) (3.16) 

+ y (x + x 3 + . . . + a 2k -iX 2k ~ 4 + x 2k ~ 3 + a 2k x 2k ). 

Both parameters act in a similar way: they rotate the vector field of (3.16) 
counterclockwise turning the origin into a nonrough unstable focus. 

Fix these parameters under «2fc-4 — a 2 k-4^ a 2k — & 2 k an d input the parameter 
«2fc-2 < into (3.16) getting again system (3.15). Since, by our assumption, 
this system has two limit cycles surrounding the origin for a 2 fc_2 > a 2 k-2> 
there exists some value of the parameter, all_ 2 (a 2k _ 2 < all_ 2 < 0), for 
which a semi-stable limit cycle, Ti2, appears in system (3.15) and then splits 
into a stable cycle, r\, and an unstable cycle, T 2 , under further decreasing 
«2A;-2- The formed domain D 2 bounded by the limit cycles I\, Y 2 and filled by 
the spirals will enlarge since, on the properties of a field rotation parameter, 
the interior unstable limit cycle will contract and the exterior stable limit 
cycle Ti will expand under decreasing a 2 k- 2 - The distance between the spirals 
of the domain D 2 will naturally increase, which will prevent the appearance 
of a semi-stable limit cycle in this domain for a 2 k- 2 < ctH_ 2 [TJ|, [T5] . 
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Thus, there are no such values of the parameters, a 2k > 0, a 2k _ 2 < 0, and 
a 2k-4 > 0) f° r which system (3.15) would have an additional semi-stable limit 
cycle surrounding the origin. Obviously, there are no other values of the para- 
meters a 2k , ®2k-2, and c%2k-i for which system (3.15) would have more than two 
limit cycles surrounding this singular point. On the same reason, additional 
semi-stable limit cycles cannot appear around the other singularities (foci or 
nodes) of (3.15). Therefore, 2 + I is the maximum number of limit cycles in 
system (3.15). 

Suppose that system (3.15) has two limit cycles, Ft and T 2 , surrounding the 
origin and I limit cycles surrounding the other antisaddles of (3.15) (this is 
always possible if a 2 k 3> — a 2k _ 2 «2fc-4 > 0). Fix the parameters a 2 k, «2fc-2, 
a 2 k~4 and consider a more general system inputting the fourth parameter, 
a 2k - 6 < 0, into (3.15): 

x = y, 

y = -x (1 + ft x ± x 2 + h x 3 ± x 4 + . . . + f3 2l ^ x 21 - 1 ± x 21 ) (3.17) 
+ y(x + x 3 + ... + a 2k - 6 x 2k ~ 6 + x 2k ~ 5 + . . . + a 2k x 2k ). 

For decreasing a 2k ^§, the vector field of (3.17) will be rotated clockwise and 
the focus at the origin will immediately change the character of its stability 
generating a third (stable) limit cycle, T 3 . With further decreasing a 2 k-e, T 3 
will join with T 2 forming a semi-stable limit cycle, T 2 3, which will disappear in 
a "trajectory concentration" surrounding the origin; the cycle Ft will expand 
disappearing on a separatrix cycle of (3.17). 

Let system (3.17) have three limit cycles surrounding the origin: Ft, T 2 , T^. 
Could an additional semi-stable limit cycle appear with decreasing a 2 k-6 after 
splitting of which system (3.17) would have five limit cycles around the origin? 
It is clear that such a limit cycle cannot appear either in the domain D 2 
bounded by the cycles Ft and T 2 or in the domain D 4 bounded by the origin 
and r 3 because of the increasing distance between the spiral coils filling these 
domains after decreasing a 2 k^§. Consider two other domains: D\ bounded on 
the inside by the cycle Ft and D 3 bounded by the cycles T 2 and T 3 . As before, 
we will prove the impossibility of the appearance of a semi-stable limit cycle 
in these domains by contradiction. 

Suppose that for some set of values of the parameters a 2k > 0, a* 2k _ 2 < 0, 
a 2k~4 > 0) an d a 2k-e < such a semi-stable cycle exists. Return to sys- 
tem (3.10) again, input first the parameters «2fc-6 < 0, a 2 t-2 < and then 
the parameter a 2k > : 
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y = -x (I + fax ± x 2 + f3 3 x 3 ± x 4 + . . . + (3 2l ^ x 21 - 1 ± x 21 ) (3.18) 
+ y (x + x 3 + . . . + a 2 k-6X 2k ' 6 + . . . + a 2 k-2X 2k ~ 2 + x 2k ~ 3 + a 2 k x 2k ). 

Fix the parameters a 2/ t_ 6 , a 2k _ 2 under the values a 2k _ 6 , & 2k - 2 , respectively. 
With increasing a 2k , a separatrix cycle formed around the origin will generate 
a stable limit cycle, r\. Fix a 2k under the value a 2k and input the parameter 
«2fc-4 > into (3.18) getting system (3.17). 

Since, by our assumption, (3.17) has three limit cycles for a 2k -4 < a 2k _ 4 , 
there exists some value of the parameter a||_ 4 (0 < a|f_4 < a 2k -4) for which 
a semi-stable limit cycle, r 2 3, appears in this system and then splits into an 
unstable cycle, r 2 , and a stable cycle, r 3 , with further increasing a 2 fc-4- The 
formed domain D 3 bounded by the limit cycles r 2 , T 3 and also the domain 
Di bounded on the inside by the limit cycle T\ will enlarge and the spirals 
filling these domains will untwist excluding a possibility of the appearance of 
a semi-stable limit cycle there [H], [T5] . 

All other combinations of the parameters a 2k , a 2k _ 2 , a 2k ^4, and tx 2k -§ are 
considered in a similar way. It follows that system (3.17) can have at most 
3 + / limit cycles. 

If we continue the procedure of successive inputting the even parameters, 
a 2k , . . . , a 2 , a , into system (3.10), it is possible first to obtain k limit cycles 
surrounding the origin (a 2k > -a 2k _ 2 > a 2k _ 4 > -a 2fe _ 6 > «2fc-8 > • • •) 
and then to conclude that the canonical system (3.5) (i.e., the Lienard poly- 
nomial system (1.3) as well) can have at most k + l limit cycles, k surrounding 
the origin and / surrounding one by one the antisaddles (foci or nodes) of (3.5) 
(and (1.3) as well). The theorem is proved. □ 
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